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Abstract : This thesis deals with the computational complexity of the 
maximum network flow problem. Ue first introduce the basic concepts and 
fundamental theorems upon which the study of "max-flow" has been built. Ue 
then trace the development of max-flow algorithms from the original 
"labeling algorithm" of Ford and Fulkerson, through a recent 0(V-E-log V) 
algorithm due to Galil and Naamad. Ue include a description of each of 
these algorithms, along with a proof of correctness and proof of running 
time for most of them. Finally we turn our attention to the problem of 
establishing lower bounds on the complexity of max-flow. Ue show that a 
straightforward application of the polyhedral lower bound technique 
developed by Yao, Avis and Rivest fails to produce a non-linear lower bound 
on max-flow. In the process, however, we prove several interesting results 
concerning the facial structure of a class of polyhedra very closely 
related to the maximum network flow problem. 
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CHAPTER 1 - INTRODUCTION 

1.1 An Overview of the Thesis 

The problem of determining a maximum steady state flow from one 
point to another in a network with edge capacity constraints, has come to 
be known as the maximum network flow problem ("max-flow" in short). L.R. 
Ford and D.R. Fulkerson [91 were the first to study max-flow as a 
computational problem. They developed the first max-flow algorithm in the 
mid 1950* s and laid the groundwork for much of the research that was to 
follow. Since that time max-flow has been widely studied and has developed 
great practical application, especially in the analysis of transportation 
and communication networks. This thesis will be concerned with the 
computational complexity of the maximum network flow problem, investigating 
both upper and lower bounds on the problem. 

The remainder of this chapter will be devoted to a development of 
the foundation necessary for any coherent study of the maximum network flow 
problem. Ue will begin by presenting the basic definitions and concepts 
that have become standard in the max-flow literature. Ue shall then 
introduce the fundamental theorems upon which the study of max-flow has 
been built. These theorems, due originally to Ford and Fulkerson [9], are 
known as the Augmenting Path Theorem, the Integral Flow Theorem and the 
Max-Flow n in-Cut Theorem. 

The second chapter in this thesis will deal with upper bounds on 
the complexity of the maximum network flow problem. The long and 
intriguing history of the search for such bounds is summarized below in 



Table 1.1. 



Upper Bounds on Max-Flow 



1) Ford and Fulkerson (1956) Unbounded 

2) Edmonds and Karp (1969) CHV-E 2 ) 

3) Dinic (1970) 0(V 2 -E) 

4) Karzanov (1973) 0(V 3 ) 

5) Cherkasky (1976) 0(V 2 -E 1/2 ) 
G) Gali I (1978) 0(V 5/3 E 2/3 ) 

7) Malhotra, Kumar and Maheshwari (1978) 0(V 3 ) 

8) Gali I and Naamad (1978) 0(V£-log 2 V) 

Table 1.1 



Each of these bounds has been demonstrated by the construction of a max- 
flow algorithm which has the specified worst case running time. Chapter 2 
will consist of a survey of each of these algorithms. Ue note here that 
the algorithm developed by Malhotra, Kumar and Maheshwari [16], which we 
will call the MKM algorithm, does not result in an asymptotic improvement 
over the previous three algorithms. However, this algorithm is very simple 
and is probably the best algorithm to use on dense networks (E«V ) . 

The determination of lower bounds on the computational complexity 
of the maximum network flow problem has thus far received little attention 
in the literature. In Chapter 3, however, we shall investigate one 
particular approach to establishing a non- linear lower bound on the 
complexity of max-flow. The technique we shall deal with is the polyhedral 
technique developed in 1977 by A.C. Yao, D.M. Avis and R.L. Rivest [211. 
Ue will show that a straightforward application of this technique fails to 



produce a non-linear lower bound on max-flow. In the process, however, we 
shall answer several questions concerning the facial structure of a class 
of polyhedra very closely related to the max-flow problem. 



1.2 Basic Definitions and Concepts 

This section will be devoted to a presentation of the basic 
definitions and concepts that have become standard in the study of max- 
flow. He will begin by defining a network and a legal flow function on a 
network. These concepts will then be used to develop a. formal definition 
of the maximum network flow problem. Finally we will introduce the basic 
notions of a cut and a flow augmenting path. 

De f i n i t i on 1.1 ; 

A network 71 = (G,s,t,c) is a 4-tuple with the properties; 

1) G = (V,E) is a finite directed graph composed of a 
set of vertices, V, and a set of edges, E. 

2) Two distinct vertices seV and teV are specified as 
the source and sink respectively. 

3) Each edge eeE is assigned a non-negative real number 
c (e) , ca I I ed the capaci ty of edge e. 

For any vertex veV, In(v) denotes the set of all edges incoming to 
v and Out(v) denotes the set of all edges outgoing from v. An edge e, 
directed from a vertex u to a vertex v in a network Tl, will often be 
denoted by the ordered pair e«=(u,v). 



Example 1.1 : (Network Rq) 

V = |s,a,b,c,d, t} 

E = {(s,a) , (s,b) , (a,c) , (a,d) , 

(b,d), (c,t), (d,t)} 
c(e) = 18 [for all eeE s.t. e*(a,d)] 
c(a,d) - 1 

In(t) - {(c,t), (d,t)} 
Out(s) = {(s,a), (s.b)} 




Def ini t i on 1.2 ; 

A legal f Ion function f on a network 71 associates with each edge 

eeE a real number f(e) satisfying the conditions: 

CI) 0<f(e)<c(e), for each edge eeE 

C2) £ He) - £ f(e) - » for each vertex veV-{s,t}. 
eeln(v) eeOut(v) 

The value of a legal flow function f is defined to be: 

(1.1) v(f) = J} f(e) ~ £ f( e). 

eeln(t) eeOut(t) 

Informally we say f is a steady state flow from s to t and f(e) is 
the steady state flow through edge e. Thu9, condition CI tells us that the 
steady state flow through any edge must be non-negative and not exceed the 
capacity of that edge. Further, condition C2 states that any steady state 
flow from s to t must have the property that flow is conserved at every 
vertex other than the source and the sink. 



Def ini t ion 1.3 : 

An edge e, incident upon vertices u and v in a network 71, is said 
to be useful from vertex u to vertex v with respect to a legal flow 



function f defined on Tl, if either: 



1) e«(u,v) and f(e)<c(e) 

or 

2) e«(v,u) and f(e)>0. 



Example 1.2 ; 

The following legal flow function f is defined on network Tin from 
Example 1.1; 

18,4 



fr»(s,a) 



r 8 



5 f (s,b) 

ig(a,c) = 4 f (a,d) 

f (b,d) = 8 f (c,t) 

f (d,t) - 1 



8 
1 
4 




18,8 



v(f,,) = fo(c,t) + fp(d,t) - 8 - 4+1 - 5 



'8 



8' 



Note 



(s,b) is useful from s to b but not useful from b to s 
(a,d) is useful from d to a but not useful from a to d 
(a,c) is useful from a to c and from c to a 



Let us now consider a legal flow function, defined on a network Tl, 
whose value is maximum over the set of values of all legal flow functions 
defined on Tl. Such a flow function is said to be maximum with respect to 
Tl. Notice that the existence of a maximum flow function on a network is 
not open to question since a network is composed of only finite capacity 
edges. Further, notice that there may exist more than one distinct maximum 
flow function on a network. 



De f i n i t i on 1 . 4 ; (flax-Flow) 

The maximum network f low problem is defined as the problem of 
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computing a maximum f low function on a network given as input. 

He now turn our attention to the basic concepts underlying the 
theory on which the study of the maximum network flow problem has been 
built. In particular, we present the notions of a cut and a flow 
augmenting path. 

Def ini t i on 1.5 : 

A cut in a network 77 is a set of vertices X with the properties: 

1) XcV. 

2) seX. 

3) teX" where R - V-X. 

The set of all edges e€E which are directed from a vertex in X to a vertex 

in R is denoted by (X;X). The capaci ty of a cut X is defined to be: 

(1.2) C(X) = £ c(e). 

ee(X;X) 



Example 1.3 ; 

The following cut Xo is defined on network 77g from Example 1.1: 



{s.a.d} 



R = {b.c.t} 




C(X ) = c(s,b) + c(a,c) + c(d,t) - 18+10+10 « 30 
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Defihi t ion 1.6 ; 

A path p, from vertex Vi to vertex v n in a network Tl, is a 
sequence of distinct vertices and edges p - v l e l v 2 e 2 "*' v n-l e n-l v n ^ n -^ 
such that VjeV (for each i-l,...,n), BjeE (for each i-1, . . . ,n-l) and 
either e ( = (v- .Vj + ^), in which case ej is said to be a forward edge in p, 
or ej = (vj + i.Vj)t in which case e- ( is said to be a reverse edge in p, for 
each i=l,...,n-l. Any path composed of only forward edges is called a 
chain . Note that a path is defined to be acyclic. 

Def ini t ion 1.7 ; 

A f low augmenting path , with respect to a legal flow function f on 
a network Tl, is a path p' with the property that each edge ej€p' is useful 
from vertex v,ep' to vertex Vj + i€p'. A flow augmenting path from the 
source to the sink is called an s-t flow augmenting path. We shall see in 
the next section that an s-t flow augmenting path is a path along which f 
can be augmented (i.e. the value of f can be increased). 

Example 1.4 ; 

The following s-t flow augmenting path pg is defined with respect 
to the legal flow function fo from Example 1.2; 

p0 = s(s,b)b(b,d)d(a,d)a(a,c)c(c, t) t 

Note - (s,b) , (b,d), (a,c) and (c,t) are all forward edges along p 
(a,d) is a reverse edge along pg. 
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1.3 Fundamental Theorems 

The study of the maximum network flow problem has been based on 
three fundamental theorems first proven by Ford and Fulkerson [9] in the 
mid 1950' s. These theorems are known as the Augmenting Path Theorem, the 
Integral Flow Theorem and the Max-Flow Hin-Cut Theorem. This section will 
be used to develop formal proofs of each of these theorems. Ue will begin 
by proving the following principle lemma which we alluded to at the end of 
the last section. 

Lemma 1.1 ; 

If there exists an s-t flow augmenting path p', with respect to a 
legal flow function f on a network Tl, then there exists a legal flow 
function f' on R such that v(f')>v(f). 

Proof ; 

Let P' =v i e i v 2 e 2' ' ,v n-l e n-l v n * n *2) be an s-t flow augmenting path 
with respect to a legal flow function f on a network Tl. Define the sets 
Eg, E^ and En as follows: 

Eg = {eeE | e is not an edge in p'} 

Ei = {eeE | e is a forward edge in p'} 

E2 » {eeE I e is a reverse edge in p'}. 
Note that E^U^ is the set of all edges in p' and that EgUEjl^ is the set 
of all edges in 77 (i.e. EgUEjl^ - E). 

Since p' is an s-t flow augmenting path, it follows from Definition 
1.7 that vj = s, v n -t and ej is useful from v- to Vj + ^, for each i-l,...,n-l. 
More explicitly, applying Definition 1.3 yields: 
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(1.31 f(e)<c(e) =» c(e)-f(e)>8, for each edge eeEi 

and 

(1.4) f(e)>8, for each edge eeEo. 
Now construct the function f as follows: 

(1.5) c(e)«-c(e)-f (e) , for each eeEi 

(l.G) c(e)<-f(e), for each eeE? 

o ♦- min c(e) 
eeEjl^ 

(1.7) f'(e)«-f(e), for each eeEg 

(1.8) f'(e)«-f (e)+i, for each eeEj 

(1.9) f'(e)«-f (e)-i, for each eeE 2 . 

It must now be shown that f is a legal flow function on Tl and that 
v(f')>v(f). (Notice that f' is defined on all edges eeE.) 

To show that f' is a legal flow function on Tl, we must simply 
prove that f' satisfies both condition CI and condition C2 of Definition 
1.2. Let us first consider condition CI. Since f is a legal flow function 
on Tl, it follows from Definition 1.2 that: 

(1.10) 0<f(e)<c(e), for all eeE. 
Thus, applying (1.7) yields: 

(1.11) 0<f'(e)<c(e), for all eeE . 

He now notice, from (1.3) and (1.4), that c(e)>8 for all eeEjl^ and thus 
fl>0. Combining this with (1.8) and (1.9), we have: 

(1.12) f'(e)>f(e), for all eeEj 

and 

(1.13) f'(e)<f(e), for all eeE 2 . 
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lie can now apply (1.10) and (1.12) to (1.8) to obtains 

(1.14) 0<f'(e)<c(e), for all eeEp 
Similarly, we can apply (1.18) and (1.13) to (1.9) to obtains 

(1.15) 0<f' (e)<c(e), for all eeE 2 . 
Therefore, combining (1.11), (1.14) and (1.15), we have: 

0<f'(e)<c(e), for all eeE 
and thus f' satisfies condition CI of Definition 1.2. 

Ue must now prove that f' satisfies condition C2 of Definition 1.2. 
This is accomplished by considering any vertex v'€V-{s, t}. Since f is a 
legal flow function on 77, it follows from Definition 1.2 that: 

(1.16) 2 f( e) - 2 f(e) " 0- 

eeln(v') eeOut(v') 

Further, if v' is not a vertex along the path p', then all of the edges 

incident upon v' must be contained in the set Eg. Thus, combining (1.7) and 

(1.16) yields: 

2 f'(e) - 2 f ' {e) - 2 f(e) ~ 2 f{e) " - 
eeln(v') eeOut(v') eeln(v') eeOut(v') 

(for all v'eV-{s,t} and v'*p') 

If, however, v' is a vertex along the path p', then let i be the index of 
v' along p' (i.e. v' = v ( ep'). Since p' is acyclic and v'^{s, t}, there 
must be exactly two distinct edges incident upon v' which are contained in 
p", namely ej_^ and ej. All other edges incident upon v' are therefore 
contained in the set Eg. He must now consider each of the following four ' 
possible cases: 
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1) ej^eEj and e^eE^ 

2) ej_^eE^ and B\^.j 

3) ej_^eE2 and ejeE^ 

4) ejieEo and BjGEo. 

For case 1, we combine (1.7), (1.8) and (1.16) to obtain: 

E f'<e> - E f'(e) - [(E '<•>) + «] " [(S *«•>) + »] 

eeln(v') eeOut(v') eeln(v') eeOut(v') 

- (E f(e)) - (E f(e)) + b - b 
eeln(v') eeOut(v') 

- E f(e) " E f{e) - - 
eeln(v') eeOut(v') 

The same result is proven for cases 2, 3 and 4 similarly. He therefore 

have that: 

E f'(e) - E f ' (e) - . for all v'eV-{s,t} 
eeln(v') eeOut(v') 

and thus f' satisfies condition C2 of Definition 1.2. 

To show that v(f')>v(f), we first notice that there is exactly one 

edge incident upon t which is contained in the path p', namely e n _j. All 

other edges incident upon t are therefore contained in the set Eo. He must 

now consider each of the following two possible cases: 

1) VieE! 

2) e n _ ie E 2 . 
For case 1, we can simply combine (1.1) of Definition 1.2, (1.7) and (1.8) 
to obtain: 

v(f') - E f ' (e) ~ E f ' (e) 
e€ln(t) ecOut(t) 



[(E f(e)) + l] - E "el 

eeln(t) eeOut(t) 
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- S f(e) - Z f(e) + * 

e€ln(t) eeOut(t) 

- v(f) + 6 > v(f). (Bince 6>0) 

The same result is proven for case 2 similarly. Ue therefore have that 
v(f') > v(f). 

D 

The proof of Lemma 1.1 reveals a fairly simple procedure for 
augmenting an existing legal flow function on a network, given a 
corresponding s-t flou augmenting path. The procedure consists of 
determining the "excess capacity" along each edge in the flow augmenting 
path, as defined by (1.5) and (1.6), and then increasing the flow along all 
forward edges and decreasing the flow along all reverse edges in the path. 
The amount by which flow is increased or decreased along each edge is 
simply the minimum excess capacity over all edges in the flow augmenting 
path. 

We now present a lemma which will be useful in demonstrating a 
relationship between the capacity of a cut and the value of a legal flow 
function on a network. 

Lemma 1.2 ; 

Given any cut X and any legal flow function f on a network Tl, 

v(f) - 2 f(e) " X f(e) - 
ee(X;5<) ee(S<;X) 

Proof ; 

Let X be any cut and f be any legal flow function on a network 77. 
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Then by condition C2 of Definition 1.2 and the fact that s $ R, we have: 

(1.17) 2 (Z f(e) " 2 f(e) ) " - 
veX-m eeln(v) eeOut(v) 

Combining (1.1) of Definition 1.2 and (1.17) yields: 

(1.18) v(f) . £ f(e) " 2 f<e) + 2 (2 f(e) " 2 Me)). 

eeln(t) eeOut(t) veX-lt) e€ln(v) eeOut(v) 

Simplifying (1.18) we obtain: 

v(f) = 2 (D f(e) - 2 f(e)) 

veR eeln(v) eeOut(v) 

- 2 2 f(e) - 2 2 f(e) - 

veR eeln(v) veR eeOut(v) 

Thus, we have that the value of f is equal to the sum of the flow along all 

edges directed into a vertex in X minus the sum of the flow along all edges 

directed out of a vertex in R. If we now consider separately those edges 

directed from a vertex in X to a vertex in R and those edges directed from 

a vertex in X to a vertex in X, we obtain: 

v(f) = (2 Me) + 2 f(e)) - (2 Me) + 2 Me)) 
ee(X;X) ee(R;R) ee(R;X) ee(R;R) 

=■ S f(e) ~ 2 f(e) - 
ee(X;R) ee(R;X) 



Corol lary 1.1 : 

If the value of a legal flow function f' is equal to the capacity 
of some cut X' on a network 77, then f' is maximum on 77 and X' has minimum 
capacity over all cuts on 77. 

Proof : 

Let X be any cut and f be any legal flow function on the network 
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Tl. Then by Lemma 1.2 we have that: 

v(f) - £ f(e) - £ f{e, « 
ee(XjR) ee(R;X) 

Applying condition CI of Definition 1.2 we obtain: 

v(f) i 2 c(e) 
ee(X;5<) 

and thus by (1.2) of Definition 1.5: 

(1.19) v(f) < C(X). 

Therefore, if (1.19) holds by equality for some legal flow function f' and 

some cut X' on the network Tl, then f' mu9t be maximum on Tl and X' mu9t 

have minimum capacity over all cuts on Tl. 

D 

14e are now ready to present the three fundamental theorems upon 
which the study of max-flow has been built. 

Theorem 1.1 : (Augmenting Path Theorem) 

A legal flow function f on a network Tl is maximum if and only if 
there exists no s-t flow augmenting path with respect to f on Tl. 

Proof : 

Clearly, if there exists an s-t flow augmenting path with respect 
to a legal flow function f on a network Tl, then by Lemma 1.1 f is not 
maximum on Tl. Assume now that there is no s-t flow augmenting path with 
respect to f on Tl and define the set S as follows: 

S = {veV | 3 an s-v flow augmenting path with respect to f on Tl\ U {s}. 
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Since there is no s-t flow augmenting path with respect to f on 77, it can 
easily be seen from Definition 1.5 that S forms a cut in 77. Further, from 
Definition 1.3 and Definition 1.7 we have that f(e)«c(e) for each edge 
ee(S;S) and f(e) =0 for each edge ee(5;S). He can now apply Lemma 1.2 to 
obtain: 

v(f) - ]£ f(e) " X f(e) 
ee(S;3") ee(§;S) 

- £ c(e) 
ee(S;S) 

- C(S). 

Thus, by Corollary 1.2 we have that f is maximum on 77. 



Theorem 1.2 : (Integral Flow Theorem) 

There exists an integral valued maximum flow function on any 
network defined by an integral valued capacity function. 

Proof : 

Let 77 be any network defined by an integral valued capacity 
function and let fx be the zero flow function on 17, defined by f^teJ-B for 
each edge eeE (Notice that such a flow function will be a legal flow 
function on any network). He can now compute a maximum flow function on 77 
as fol lows: 

wh i I e there exists an s-t flow augmenting path 
with respect to fa on 77 

do augment fa as outlined in 
the proof of Lemma 1.1. 
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An examination of (1.5) through (1.9) reveals that the legal flow function 
generated at each iteration of this procedure will be integral valued. 
Thus by Lemma 1.1, the value of each successive legal flow function 
generated must be at least one integral unit greater than the value of the 
previous legal flow function. Combining this fact with (1.19) and Theorem 
1.1, we now have that our procedure must halt within a finite number of 
steps, yielding an integral valued maximum flow function on Tl. 

D 

Theorem 1.3 ; (Max-Flow M in-Cut Theorem) 

The value of any maximum flow function on a network Tl is equal to 
the minimum cut capacity over all cuts on Tl. 

Proof ; 

Let f be any maximum flow function defined on a network Tl. 
Applying Theorem 1.1 we have that there is no s-t flow augmenting path with 
respect to f on Tl. 1 1 now follows immediately from the proof of Theorem 
1.1 that there exists a cut S on Tl such that: 

v(f) - C(S). 
Further, by Corollary 1.2 we have that the cut S must have minimum capacity 
over all cuts on Tl. Thus, the value of any maximum flow function on Tl is 
equal to the minimum cut capacity over all cuts on Tl. 

D 

Ue shall see in Chapter 2 how the previous three theorems form the 
basis for all the max-flow algorithms that have thus far been developed. 
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In fact, we will see that each algorithm is actually a variation of the 
procedure given in the proof of Theorem 1.2. 
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CHAPTER 2 - UPPER BOUNDS ON MAX-FLOW 

2.1 Introduction 

The formal definition of the maximum network flow problem, as 
presented in Section 1.2, can be explicitly stated as follows. Given any 
network Tl as input, compute values of the variables X_ [for each eeE] so 
as to maximize the objective function 

S x e - 2 x e 

eeln(t) eeOut(t) 



subject to the constraints 



X e > 0, for each eeE 

X e < c(e), for each eeE 

£ x e - £ x e = 6, for each veV-{s,t}. 
eeln(v) eeOut(v) 

Thus the maximum network flow problem can be viewed as an optimization 

problem in which a linear function must be maximized subject to a system of 

linear equations and linear inequalities. G.B. Dantzig [4] developed an 

algorithm in the early 1950* s, known as the simplex method, which could be 

used to solve such linear programming problems. Although it would not be 

incorrect to consider the simplex method to be the first max-flow 

algorithm, it is usually not treated as such. The simplex method is a very 

general a I gor i thm which has an unbounded worst case running time. He have 

included it here only for the sake of completeness. 

L.R. Ford and D.R. Fulkerson [9] were the first to produce 

significant research results concerning, specifically, the maximum network 

flow problem. In the mid 1950's they proved the Augmenting Path Theorem, 
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the Integral Flow Theorem, and the flax-Flow Min-Cut Theorem. These 
fundamental results led directly to their development of the labeling 
algorithm for solving the maximum network flow problem. The labeling 
algorithm is a straightforward algorithm which simply augments an existing 
legal flow function along some s-t flow augmenting path in a network. This 
process is then repeated until there no longer exist any s-t flow 
augmenting paths in the network. The labeling algorithm, although it also 
has an unbounded worst case running time, remained in successful use for 
almost 15 years. 

In 1969, J. Edmonds and R.N. Karp t7] developed a variation of the 
labeling algorithm which utilized a Breadth First Search in picking out the 
s-t flow augmenting paths in order of increasing length. This resulted in 
a much more efficient algorithm with a bounded 0(|V|-|E| ) worst case 
running time. Independently and a short time later, E.A. Dinic [5] 
developed an improved version of Edmonds and Karp's algorithm. Dinic also 
utilized the technique of Breadth First Search but he developed an 
.algorithm with time complexity 0(|V| Z -|E|). 

A.V. Karzanov [14] modified Oinic's algorithm in 1973 to obtain an 
( | V | ) max- flow algorithm. Karzanov' s algorithm was unique in that it 
simultaneously augmented an existing legal flow function along several s-t 
flow augmenting paths. In 197G, B.V. Cherkasky [183 showed how to combine 
Oinic's algorithm with Karzanov's algorithm to produce a new and very 
complex 0(|V| 2 -|E| 1/2 ) max-f low algori thm. Two years later, Zvi Galil [10] 
improved Cherkasky' s algorithm to ( | V | 5/r3 - |E | 2/3 ) by developing a 
technique for retaining useful information about the structure of the 
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network. 

In 1978 V.M. Malhotra, fl. Pramodh Kumar and S.N. Maheshwari tlG] 
discovered a very simple 0(|V| 3 ) max-flow algorithm which we shall call the 
MKtt algorithm. Their algorithm, similar to Karzanov's algorithm in that it 
simultaneously augments along several s-t flow augmenting paths, requires 
very little overhead. Although the flKfl algorithm does not result in an 
asymptotic improvement over the previous three algorithms, its simplicity 
makes it perhaps the best algorithm to use on very dense networks (E*V Z ) . 

Finally, Gal i I and A. Naamad [12] have recently developed a 
modification to the original Dinic algorithm which results in an algorithm 
with time complexity ( | V | -|E| -log 2 | V | ) . Their modification of Dinic's 
algorithm is similar to Galil's modification of Cherkasky's algorithm. 
Once again a technique is developed for retaining useful information about 
the structure of the network. 

The remainder of this chapter will be devoted to a closer 
examination of each of the maximum network flow algorithms. 

2.2 Ford and Fulkerson 

Ford and Fulkerson [9] developed the first maximum network flow 
algorithm, known as the labeling algorithm, in 195G. Their algorithm simply 
augments, along some s-t flow augmenting path, the existing legal flow 
function on a network. This augmentation is then repeated until there no 
longer exist any s-t flow augmenting paths on the network. In practice the 
zero flow function (i.e. f(e)-0 for all edges e) is used as the initial 
existing legal flow function. 
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The labeling algorithm i9 composed of two basic routines which are 
iterated until a maximum flow function is computed. The first routine 
essentially searches in a systematic way for an s-t flow augmenting path on 
the network. The second routine then augments the existing legal flow 
function along this flow augmenting path. The actual flow augmentation is 
performed exactly as outlined in the proof of Lemma 1.1. The following 
explanation of the labeling algorithm is taken from [133 : 

Step 1. Labeling Process. 

Every vertex is always in one of three states, labeled and scanned, 
labeled and unscanned, or unlabeled. A vertex is labeled and scanned if it 
has a label and we have inspected all vertices adjacent to it. A vertex is 
labeled and unscanned if it has a label but not all vertices adjacent to it 
have been inspected. A vertex is unlabeled if it ha9 no label. 

Initially, all vertices are unlabeled. A label for a vertex v. 
always has two parts. The first part is the index of a vertex Vj, which 
indicates that we can send flow from v- ( to v., and the second part is a 
number which indicates the maximum amount of flow we can send from the 
source to v. without violating the capacity constraints. He first assign 
the label [s + ,€(s)=eo] to the source, v_. The first label simply says that 
we can send flow from the source to itself; the number » indicates that 
there is no upper bound on how much can be sent. The source is now labeled 
and unscanned and all other vertices are unlabeled. In general, select a 
vertex v. which is labeled and unscanned. Assume v- has a label of the 
form ti + ,e(j)] or ti~,e(j)]. For all adjacent vertices v^ which are 

unlabeled, adjacent to v- via an edge directed from v- to v u , and for which 

J J * 



26 

the edge e=(v., v k ) is useful from v. to v^ (i.e. f (v-, v k )<c(v.,v k ) ) , assign 
the label [j + ,e(k)I to v k , where: 

e(k) =min[e(j), c(v., v k )-f (v-,v k )] . 
For all adjacent vertices v k which are unlabeled, adjacent to v- via an 
edge directed from v k to v-, and for which the edge e=(v k ,v.) is useful 
from v. to v k (i.e. f(v k ,v.)>0), assign the label [j~,e(k)l to v k , where: 

elk) - minte(j), f(v k ,v-)]. 

The + and the - signs in the labels indicate whether the 
corresponding edges appear as forward or reverse edges in the s-t flow 
augmenting path. Now all the vertices adjacent to v. have labels; v. is 
considered to be labeled and scanned and may be disregarded during the rest 
of this step. (If one inspects all the vertices adjacent to Vj and cannot 
label all these vertices, then v. is also considered to be a labeled and 
scanned vertex.) All the vertices v k are now labeled and unscanned. 

Continue to assign labels to vertices adjacent to labeled and 
unscanned vertices until either the sink is labeled or no more labels can 
be assigned and the sink is unlabeled. If the sink cannot be labeled, no 
s-t flow augmenting path exists and, hence, the existing flow function is 
maximum. If the sink is labeled, an s-t flow augmenting path has been 
found and the flow augmentation can be performed using step 2. 

Step 2. Flow Change. 

Assume that the sink is labeled tk + ,e(t)]. Let 

f (v k ,v t )«-f (v k ,v t )+e(t) and turn to v R . If v k is labeled [j + ,€(k)l, let 

f (v.,v k )«-f (v.,v k )+e(t) and turn to v.. If v k is labeled [j~,e(k)], let 

f (v. , v .)«-f (v., v -)-e(t) and turn to v.. Continue until the source is 
K J * J J 
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reached. Erase the labels on all the vertices and go back to step 1. 

Uhen the labeling algorithm terminates, the set of all labeled 
vertices clearly forms a cut in the network. Further, applying Lemma 1.2 
reveals that the capacity of this cut must be equal to the value of the 
existing legal flow function on the network. Thus by Corollary 1.1, the 
existing legal flow function must be maximum. It now remains to be shown 
that the labeling algorithm will always terminate within a finite number of 
steps. 

Let us first consider the case in which the network is defined by 

an integral valued capacity function. By the same argument as that used in 

the proof of Theorem 1.2, it can easily be seen that the labeling algorithm 

will terminate after at most v(f M . w ) iterations, where v(f M _„) is the 

max max 

f i n i t e value of a maximum flow function on the network. Further, each 
iteration will require at most 0(|E|) operations since each edge is 
examined at most twice in the labeling procedure and at most once in the 
augmenting procedure. Thus we have that the labeling algorithm will 
correctly compute, in time ( |E | -v tf ma>< ) ) , a maximum flow function on any 
network defined by an integral valued capacity function. It can also be 
shown, however, that there actually exist networks which force the labeling 
algorithm to perform v(f m ) iterations. Consider for example the network 
H, in Figure 2.1. 



Network Ri 
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Figure 2.1 



If the labeling algorithm, beginning with the zero flow function on 77^, 
augments only along the paths s (s, a) a (a, c) c (c, f ) f (f , h) h (h, t) t and 
s(s,b)b(b,d)d(d, f) f (c, f)c(c,e)e(e,g)g(g, t) t in alternating order, then 2-N 
" v{f max* flow augmentations will be required. Thus the algorithm will 
iterate v(f m ) times. Notice that the inefficiency in this example is 
based on the fact that the labeling algorithm permits the augmentation, at 
each iteration, along any one of several existing s-t flow augmenting 
paths. 

Ue shall now consider the case in which the network is not defined 
by an integral valued capacity function. Ford and Fulkerson [9] were able 
to demonstrate the somewhat surprising result that their labeling algorithm 
might fail to terminate if the network was composed of irrational edge 
capacities. This result was based on interpreting the labeling process 
broadly enough to permit the selection of any s-t flow augmenting path at 
each iteration of the computation. Thus the labeling algorithm essentially 
has an unbounded worst case running time. Ue remark, however, that 
computers only deal with rational numbers and thus in practice we could 
expect the labeling algorithm to halt and yield a correct answer. In fact, 
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despite its weaknesses, the labeling algorithm was successfully used for 
almost IB years. 

2.3 Edmonds and Karp 

In 19B9 Edmonds and Karp [7] showed how the labeling algorithm 
could be modified to obtain a bounded worst case running time. In light of 
our previous remarks, it should not be surprising to learn that their 
modification was essentially an ordering on the selection of the s-t flow 
augmenting paths. Edmonds and Karp suggested augmenting along the shortest 
s-t flow augmenting path (i.e. an s-t flow augmenting path containing a 
minimum number of edges) at each iteration. This can be easily 
accomplished by modifying the labeling process so that the vertices are 
scanned in the same order in which they receive labels (i.e. by imposing a 
Breadth First Search on the labeling process). The remainder of the 
labeling algorithm is unchanged. The running time bound on Edmonds and 
Karp's "first labeled, first scanned" modification of the labeling 
algorithm is derived from the following results [15]. 

Consider any network R upon which there is defined a legal flow 

(k) 
function f. Let <r u denote the minimum number of edges in an s-u flow 

(k) 

augmenting path after k augmentations of f. Similarly, let r y denote the 
minimum number of edges in a u-t flow augmenting path after k augmentations 
of f. 



Lemma 2.1: 



If each flow augmentation of f is made along an s-t augmenting path 
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with a minimum number of edges, then: 

^(k+l) ^ „(k) 
u u 
and 

-(k+1) > -(k) 



for al I u,k. 

Proof ; (From [151) 

Assume that ff u < * u k , for some u,k. Moreover, let: 

(2.1) (f u k+1) = min[^ k+1) | ,]*+" <*< k) }. 

Clearly <r (k+1) > 1 (only ffl k+1) - 0), and there must be some final edge 

» u 3 s 

(u,v) or (v,u) in a shortest s-u flow augmenting path after the (k+1) 
augmentation of f. Suppose this edge is (v,u), a forward edge, with 
f (v,u)<c(v,u) (the proof is similar for (u,v)). Then o^ k+1) - Cy k+1 ' + 1 
and by (2.1) , 

(2.2) (r u k+1) * <rj k) + 1. 

Further, it must have been that f (v,u)«c(v,u) after the k th augmentation of 
f; otherwise <r u k) < tr^ k) +l < <r u k+1) , contrary to the assumption. But if 
f (v,u)=c(v,u) after the k tn augmentation of f and f (v,u)<c(v,u) after the 
(k+l) st augmentation of f, it follows that (v,u) was a reverse edge in the 
(k+l) st s-t flow augmenting path along which f was augmented. Since that 

path contained a minimum number of edges, 

„(k) _ _(k) . , 
'v % + l ' 

Combining this with (2.2), however, we obtain: 

*u +2 S % 

(k+1) (k) 

contrary to our assumption. The assumption that 9 < <r u is 
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therefore false. 

The proof that r^ k+1) 2: r u k) parallels the above. 



Theorem 2.1 : 

If each flow augmentation of f is made along an s-t augmenting path 
with a minimum number of edges, then a maximum flow function is obtained 
after no more than |V|-|E|/2 augmentations of f. 

Proof : (From [151 ) 

Each time an augmentation of f is made, at least one edge in the s- 

t augmenting path is "critical" in the sense that it limits the amount of 

augmentation. The flow through such an edge (u,v) is either increased to 

capacity or decreased to zero. Suppose (u,v) is a critical edge in the 

(k+1) s-t augmenting path. The number of edges in the augmenting path is 
«r<k> + r (k) _ (k) (k) 

U U V ^ 7 v* 

The next time edge (u,v) appears in an s-t augmenting path, say the 
(£+l) s , it will be with the opposite orientation. That is, if it was a 
forward edge in the (k+l) st , it is a reverse edge in the U+l) st , and vice 
versa. If (u,v) was a forward edge in the (k+l) st s-t augmenting path (the 
proof is similar for a reverse edge), then: 



and 



«r (k) - <r (k) + 1 
v u 

% * *v + 1 - 



By Lemma 2.1, however, <r^ } > <r' k) and r}P > T, ( , k) . Thus we have that: 
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It follows that each succeeding s-t augmenting path in which (u,v) is a 
critical edge is at least two edges longer than the preceding one. 

No flow augmenting path may contain more than |V|-1 edges. 
Therefore, no edge may be a critical edge more than |V|/2 times. But each 
s-t augmenting path along which f is augmented contains a critical edge. 
Therefore there can be no more than |V|-|E|/2 successive s-t flou 
augmenting paths and this completes the proof. 

□ 

Since Edmonds and Karp's algorithm differs from the original 
labeling algorithm only in the order in which the unscanned vertices are 
scanned, it follows that their algorithm will yield a maximum flow function 
if it halts and that their algorithm will require only 0(|E|) operations 
per iteration. By Theorem 2.1, however, Edmonds and Karp's algorithm is 
guaranteed to halt after at most |V|-|E|/2 iterations. Further, this 
result is independent of the capacity function; holding for irrational 
valued capacity functions as well as integral valued capacity functions. 
Thus Edmonds and Karp's algorithm will correctly compute, in time 
0(|V|-|E| ), a maximum flow function on any network given as input. 

2.4 D i n i c 

In 1970, E.A. Dinic [5] developed a maximum network flow algorithm 
with a bounded 0(|Vp-|E|) worst case running time. Dinic's algorithm, 
like Edmonds and Karp's algorithm, is based on successive flow 
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augmentations along s-t flow augmenting paths of minimum length. Dinic, 
however, noticed that a single breadth first search of a flow network could 
be used to isolate all minimum length s-t flow augmenting paths. Based on 
this observation, Dinic's max-flow algorithm is much more efficient than 
the "first labeled, first scanned" algorithm of Edmonds and Karp. In order 
to present Dinic's algorithm, we must first introduce the notion of a 
layered network. 



Def ini t ion 2.1 : 

A layered network LN i s a network whose vertex set is partitioned 
into disjoint subsets VoiVi,...,V| such that Vo - {s} and V^ - {t}. Ue say 
that Vj is the i th layer in LN (for 0<i<£) and that £ is the length of the 
layered network LN. Each edge e=(u,v) in LN has the property that if ueV. 
then v€Vj + ^ (i.e. every edge in LN is directed from one layer to the 
next) . 



Example 2.1 ; (Layered Network LNg) 

V ={s}, V r {a,b}, V r {c,d}, V 3 -{t} 

E={(s,a) , (s,b) , (a,c) , (a,d) , 

(b,d), (c, t), (d,t)} 



c(e)-l [for al I eeE] 
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Each iteration of Dinic's algorithm is called a phase and each 
phase is divided into two procedures. The first procedure generates a 
layered network from the original input network, in such a way that the 
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layered network isolates all existing minimum length s-t flow augmenting 
paths. The second procedure then uses this layered network to successively 
augment the existing legal flow function along these minimum length s-t 
flow augmenting paths. Ue will first describe how Dinic's algorithm 
performs each of these procedures and then we will state the entire 
algori thm. 

The first task performed during each phase of Dinic's algorithm is 
the construction of the layered network. The layers composing this network 
are created from a breadth first search of the input network. This breadth 
first search begins at the source and traverses only forward directed edges 
with flow less than capacity or backward directed edges with flow greater 
than zero (i.e. only "useful" edges). The search terminates when either 
the sink is reached or no new vertices can be visited and the sink has not 
been reached. When the sink cannot be reached, however, the existing legal 
flow function is a maximum flow function and Dinic's algorithm halts. 
Notice that performing the search in this way assures that there exists a 
flow augmenting path from the source to each vertex visited. The layers, 
Vj, are finally formed by partitioning the vertices visited according to 
the length of their path of discovery from the source. Every vertex veV. 
will then have the property that the shortest s-v flow augmenting path, 
with respect to the existing legal flow function on the input network, is 
of length i. Therefore, the length of the layered network constructed will 
be equal to the length of the shortest existing s-t flow augmenting path. 
The following procedure formalizes this construction t81 . 
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procedure LN(?7, f ) 
beg i n 

i : =8; 



wh i I e ( i : « i +1 ) > do 
begin 

T:={veV | vgV. for j<i and there exists a 

useful edge from a vertex in V.j to v}j 
i f T»^ then halt (the existing f is maximum) 
else j_f_ teT then beg i n 

l: ■ i } 

. i«-i 

end 



else V • : «T 



end 



end. 



Once the layers have been generated, the layered network's edge set 
and capacity function are constructed. The edge set is simply constructed 
from all edges in the original input network which are useful from a vertex 
in Vj to a vertex in V. |+ ^ (for all 0< 1 <j£) . Every edge in the layered 
network, however, is directed from the i tn layer to the (i+l) st layer, 
regardless of its orientation in the input network. Any edge in the 
layered network whose orientation is different in the input network is said 
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to be a reverse edge in the layered network. Similarly, any edge in the 
layered network whose orientation is the same in the input network is said 
to be a forward edge in the layered network. The capacity function c is 
then created from the excess capacity along each edge in the layered 
network as follows: 

c(e)«-c(e)-f (e) , for all forward edges in the layered network 

and 
c(e)«-f(e), for all reverse edges in the layered network, 
where f is the existing legal flow function and c is the capacity function 
on the input network. It should now be clear that every s-t chain in the 
layered network corresponds to an existing minimum length s-t flow 
augmenting path in the input network. Further, every minimum length s-t 
flow augmenting path in the input network is represented by an s-t chain in 
the layered network. Thus the layered network constructed during each 
phase of Dinic's algorithm essentially isolates all existing minimum length 
s-t flow augmenting paths. 

The second task performed during each phase of Dinic's algorithm is 
the flow augmentation. Starting at the sink in the newly constructed 
layered network, the algorithm follows edges backward to the source to find 
an s-t flow augmenting path. (Recall that every s-t chain in the layered 
network corresponds to an existing s-t flow augmenting path in the input 
network.) The existing legal flow function is then augmented along this 
path as outlined in the proof of Lemma 1.1. Notice that the excess capacity 
along each edge in the flow augmenting path is given by the layered network 
capacity function, c. After augmenting along this path, the algorithm 
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adjusts the capacity function c to reflect the new excess capacity along 
each edge in the path. It also deletes from the layered network all edges 
in the path whose excess capacity drops to zero. Finally, it deletes from 
the layered network all vertices and their incident edges which are no 
longer reachable from the source or the sink. This is accomplished by 
deleting all vertices which either have no incoming edges or have no 
outgoing edges, and continuing to delete such vertices until every vertex 
left in the layered network has at least one incoming and one outgoing 
edge. After all necessary deletions have been performed, the algorithm 
searches for another flow augmenting path and the augmentation process is 
repeated. This continues until there no longer exist any s-t chains in the 
layered network. 

Dinic's complete max-flow algorithm is: (From [20]) 



procedure DlNIC(ft): 
beg i n 

initialize existing legal flow function f, on input network Tl, to 8; 
whi le "true" do 
beg i n 

construct layered network, LN; 
for each vertex v in LN do 
beg i n 

calculate indegree (v); 
calculate outdegree (v); 

i f (indegree (v)=0) or (outdegree (v)=0) then 
add v to nul I ist 
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end ; 
uhi le t is a vertex in LN do 
beg i n 

trace back from t to 3 to find an augmenting path; 
augment f along this path; 

update c with new excess capacity along each edge in path; 
delete from LN all edges along path which now have zero 
excess capacity (i.e. c(e)=8), updating indegrees, 
outdegrees, and nullist; 
wh i I e some vertex v is on nullist do 

delete v and incident edges from LN and from nullist, 
updating indegrees, outdegrees, and nullist 
end 



end 



end. 



Recall that Dinic's algorithm terminates when the breadth first 
search performed in constructing each layered network fails to reach the 
sink. When this occurs, however, the set of all visited vertices clearly 
forms a cut in the input network. Further, applying Lemma 1.2 reveals that 
the capacity of this cut must be equal to the value of the existing legal 
flow function. Thus by Corollary 1.1, the existing legal flow function 
must be maximum. It now remains to be shown that Dinic's algorithm will 
always terminate within time 0(|Vp-|E|). This is proven as a consequence 
of the following lemma [8] which shows that the number of phases is bounded 
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by I V J . 

Let fi^ denote the length of the layered network constructed during 
the k th phase of Dinic's algorithm. 

Lemma 2.2 : 

If the (k+l) st phase of Dinic's algorithm is not the last, then 

«k + i > V 

Proof : <From [8]) 

Consider any s-t chain in the layered network constructed during 
the (k+l) st phase of Dinic's algorithm: 



e l *2 % +l 



S > V, > . . . Ve i ^ t. 



First, let us assume that all the vertices in this chain appear in 
the k th layered network. Let V. be the j th layer of the k th layered 
network. Ue claim that if v a eV b then aib. This is proven by induction on 
a. For a=0, (vg=s) the claim is obviously true. Now assume v a+ i € ^ C ' If 
c<b+l the inductive step is trivial. If, however, c>b+l then the edge e g+ j 
was not used in the k phase since it was not even in the k l layered 
network, in which only edges between adjacent layers appear. But if e g+ j 
was not used in the k* phase and is useful from v g to v g+ j in the 
beginning of the (k+l) st phase, then it was useful from v g to v g+ j in the 
beginning of the k phase. Thus, v g+ j cannot belong to V c (by procedure 
LN). Now, in particular, t=Vjj and teV^ . Therefore, \ + i^^- 
Further, equality cannot hold because then the entire s-t chain would have 
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been in the k th layered network, and if all its edges are still useful at 

the beginning of the (k+l) 9t phase then we have a contradiction to the 

termination of the k phase. 

If not all the vertices in the s-t chain appear in the k layered 
e , 
network then let Vg — ^"^ v a+l be tne f ' lrst ed 9 e such that * or 9ome b * 

v a eV b but v a+l ' s not ' n tne layered network. Thus, e a+ j was not used 
in the k phase. Since it is useful in the beginning of the (k+1) 9 
phase, houever, it was also useful in the beginning of the k phase. Thus 
the only possible reason for v < not to belong to Vjj+1 is that b+l«£^. 
Further, by the argument of the previous paragraph a>b. Therefore, a+l>£^ 
and so \ + i > \« 

a 



Corol lary 2. 1 : 

Given any network R as input, the number of phases performed by 
Oinic's algorithm must be less than or equal to |V|. 

Proof ; 

Any layered network constructed by Dinic's algorithm must contain 
no more than |V| layers. Thus by Lemma 2.2 there can be at most |V| 
phases. 

a 

Ue now notice that the time required during each phase of Dinic's 
algorithm, to construct the layered network and initialize the indegrees, 
outdegrees, and nu Mist is bounded by 0(|E|). Further, each flow 
augmentation requires time 0(|V|) and there can be at most 0(|E|) such 
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augmentations since each augmentation causes the deletion of at least one 
edge from the layered network. Finally, the total time required during 
each phase of Dinic's algorithm to delete edges, delete vertices, and 
update indegrees, outdegrees, and nullist is bounded by 0(|E|). This 
results from the fact that each edge and each vertex can be deleted from 
the layered network at most once. It should now be clear that each phase 
(iteration) of Dinic's algorithm has time complexity 0(|V|-|E|). Thus by 
Corollary 2.1, we have that Dinic's algorithm will always terminate within 

2 

time 0(|Vp-|E|). Therefore, Dinic's algorithm will correctly compute, in 
time 0(|V| -|E|), a maximum flow function on any network given as input. 

2.5 Karzanov 

A.V. Karzanov [8,14] modified Dinic's algorithm in 1973 to obtain 
an 0(|V|°) maximum network flow algorithm. Karzanov noticed that the 
layered network constructed during each phase of Dinic's algorithm could be 
used to si mul taneously augment along all existing minimum length s-t flow 
augmenting paths. He then showed how this simultaneous augmentation could 
be performed in time 0(|Vp). Karzanov'9 algorithm is the result of 
replacing Dinic's 0(|V|-|E|) successive flow augmentation procedure with 
this new OtlVp) simultaneous flow augmentation procedure. 

Karzanov' s results are based on the notion of a maximal flow 
function. A maximal flow function f, on a network Tl, is defined to be any 
legal flow function on Tl which has the property that every s-t chain in Tl 
contains at least one saturated edge (i.e. at least one edge e such that 
f(e)=c(e)). From the following example it can be seen that a maximal flow 
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function on a network need not be a maximum flow function on that network. 

Example 2.2 ; 

The following maximal flow function f is defined on the layered 
network LN from Example 2.1: 



1.8 



f (s,a)=l f (s,b)«8 

f (a,c)=8 f <a,d)=l 
f (b,d)=0 f (c,t)=B 

f (d,t)=l 




v(f )-l is not maximum. 



Karzanov noticed that the flow augmentation required during each 
phase of Dinic's algorithm could be achieved by simply augmenting the 
existing legal flow function f with any maximal flow function f on the 
current layered network. Once such a maximal flow function had been 
computed, the flow modification could be performed as follows: 

f ' (e)t-f (e)+f (e) , for all forward edges in the layered network 

and 
f ' (e)«-f (e)-f (e) , for all reverse edges in the layered network. 
It can easily be seen that performing the augmentation in this manner is 
essentially the same as simultaneously augmenting along all existing 
minimum length s-t flow augmenting paths, where each individual 
augmentation is performed exactly as outlined in the proof of Lemma 1.1. 
Karzanov used this clever augmentation technique to develop the following 
modification of Dinic's algorithm. 
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procedure KARZANOV {Tl) : 
begin 

initialize existing legal flow function f, on input network Tl, to 8; 
wh i I e "true" do 
begin 

construct layered network, LN; 

A 

compute maximal flow function f on LN; 
for each edge e in Tl do 

A 

i f e is a forward edge in LN then f (e) :-f (e)+f (e) 

a 

el se i f e is a reverse edge in LN then f (e) :«f (e)-f (e) 

end 
end. 



Notice that Karzanov's algorithm, like Dinic's algorithm, 
terminates when the breadth first search performed in constructing each 
layered network fails to reach the sink. By the same argument as that used 
for Dinic's algorithm, however, the existing legal flow function must be 
maximum when this occurs. Therefore if Karzanov's algorithm halts, then 
the existing legal flow function must be a maximum flow function on the 
input network. Next we notice that the proof of Lemma 2.2 is valid for 
Karzanov's algorithm as well as Dinic's algorithm. Thus given any network 
Tl as input, the number of phases performed by Karzanov's algorithm must be 
less than or equal to |V|. It should also be clear, however, that the time 
required to construct each layered network and perform each flow 
modification is bounded by 0(|E|). Therefore if we let t denote the time 
required to compute each maximal flow function, then Karzanov's algorithm 
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is guaranteed to halt within time 0(|V|-|E| + |V|-f). Ue will now 9how how 

2 

Karzanov's algorithm computes each maximal flow function in time C | V | ) to 

yield an 0(|V|°) max-flow algorithm. 

Karzanov's algorithm computes each maximal flow function by 
successively improving an existing illegal flow function, called a pref low 
function, on the current layered network LN [8, IB, 14) . For each vertex v 
in LN, let In'(v) denote the set of all edges incoming to v in LN and let 
Out'(v) denote the set of all edges outgoing from v in LN. A pref low 
function f, on the layered network LN, associates with each edge e in LN a 
real number f(e) satisfying the conditions: 

C3) 0<f(e)<c(e), for each edge e in LN 

C4) 2} f(e) > 2} f(e), for each vertex v*s,t in LN. 
eeln'(v) eeOut'(v) 

Thus a preflow function is simply a flow function which satisfies the 
capacity constraint but not necessarily the conservation constraint of a 
legal flow function. For each vertex v*s,t in LN, we define excess (v) to 
be the excess flow entering v: 

excess(v) = ^ f(e) - £, f(e). 
eeln'(v) eeOut'(v) 

If excess (v)>0 then v is said to be unbalanced ; otherwise v is said to be 
balanced . The source and the sink are always considered balanced. 

Throughout the execution of Karzanov's maximal flow procedure, 
every edge in LN is declared either open or closed . Initially all edges 
are declared open. As the algorithm proceeds, however, some of the edges 
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will be declared closed. Once an edge is declared closed, the flow through 
it will remain unchanged to the end of the procedure. 

Karzanov's maximal flou procedure alternates between pushing 
additional flow from unbalanced vertices and balancing the unbalanced 
vertices that are generated during these pushes. The pushing of flow is 
achieved through repeated calls to a procedure PUSH(i), with increasing i 
[101. The procedure PUSH(i) considers in turn each unbalanced vertex in 
layer v\ , attempting to push flow from it to vertices in layer V j + i« For 
each unbalanced vertex veVj, the procedure considers in turn each open 
edge in Out'(v) and sends through it the maximum possible amount of flow. 
(The two constraints that exist are the current excess of v and the amount 
of flow needed to saturate the edge.) The push from v ends when either v 
becomes balanced or every edge in Out'(v) becomes either saturated (i.e. 
f(e)=c(e)) or closed. For each vertex u in LN there is a stack (push-down 
store) on which the history of additions of incoming flow into u is 
recorded. When the flow in an edge e»(v,u) is incremented by an amount 6, 
the pair (v.fl) is added to the top of the stack for vertex u. The procedure 
PUSH ( i ) is said to be successful if flow is pushed to layer V\ + i> If 
PUSH(i) is successful then PUSH(i+l) is called, and so on. 

A procedure BALANCE ( i ) is the tool through which unbalanced 
vertices become balanced. The procedure BALANCE(i) uses the stacks to 
shift back flow from vertices in layer V- to vertices in layer V._^. It 
balances in turn each of the unbalanced vertices veV. by canceling the 
most recent additions of flow into v. Clearly the last canceled addition of 
flow into v may only be partial. After each unbalanced vertex v is 
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balanced, all the edges in In'(v) are declared closed. The procedure 
BALANCE(i) is always followed by a call to PUSH(i-l). 

Karzanov's complete maximal flow procedure is: (From [18]) 



procedure MAXIMAL (LN) : 
begin 

initialize existing preflow function f, on layered network LN, to 0; 
empty the stacks of all vertices in LNj 
i : =0; 
PLOOP: PUSH(i); 

wh i I e the previous push was successful and 'i+1<jE do 
begin 

i : = i +1 ; 
PUSH(i) 
end ; 
i f there exist unbalanced vertices in LN then 
begin 

i:=number of highest layer V.(8<j<£) containing 

unbalanced vertices; 
BALANCE ( i ) ; 
i; =i-l; 
goto PLOOP 
end ; 
for each edge e in LN do f(e);«f(e) 



end. 
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It should now be clear that every vertex in the layered network LN 
will be balanced when Karzanov's maximal flow procedure halts. Thus the 
final existing preflow function f will in fact be a legal flow function on 
LN. In order to show that it will also be a maximal flow function on LN, 
we must introduce the notion of a blocked vertex. A vertex v in LN is said 
to be blocked with respect to the existing preflow function if every v-t 
chain in LN contains at least one saturated edge. Notice that the source s 
becomes blocked after the first execution of PUSH (8), since every edge in 
Out'(s) becomes saturated. 

Lemma 2.3 ; 

If a vertex in LN becomes blocked at some point in the execution of 
Karzanov's maximal flow procedure, then it remains blocked to the end of 
the procedure. (A proof of this lemma appears in [8].) 

Lemma 2.4 : 

Every vertex in LN is balanced at most once throughout the 
execution of Karzanov's procedure. (A proof of this lemma appears in [81.) 

Ue can now see that the final existing preflow function on LN will 
be a legal flow function which has the property that every s-t chain in LN 
contains at least one saturated edge. Therefore when Karzanov's maximal 
flow procedure halts, the existing preflow function f will in fact be a 

maximal flow function on the layered network LN. It now remains to be 

2 

shown that the procedure will always halt within 0(|V| ) steps. 
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The total number of steps performed by Karzanov's maximal flow 
procedure is clearly bounded by the total number of flow additions and flow 
reductions performed. The number of flow reductions performed, however, is 
bounded by the number of flow additions performed since each vertex is 
balanced at most once and the history of flow additions in the stacks is 
used to perform the flow reductions. Thus it suffices to show that the 
number of flow additions performed by Karzanov's procedure is bounded by 
0(|V| ). Ue first notice that there can be at most one saturating flow 
addition per edge in the layered network. Since the number of edges in the 
layered network is bounded by the number of edges in the original input 
network, however, there can be at most 0(|E|) saturating flow additions. 
Next we notice that there can be at most one non-saturating flow addition 
per vertex in the layered network, between any two successive calls to 
BALANCE(i). (Uhen flow is pushed from a vertex v in layer V-, only the 
last edge considered in Out'(v) does not necessarily become saturated.) 
Since the number of vertices in the layered network is bounded by the 
number of vertices in the original input network, however, there can be at 
most ( | V | ) non-saturating flow additions between any two successive calls 
to BALANCE(i). From Lemma 2.4, however, there can be at most 0(|V|) calls 

to BALANCE(i). Thus the total number of non-saturating flow additions is 

2 
bounded by 0(|V| ) and hence the total number of flow additions is bounded 

by 0(|V| 2 ). 

It should now be clear that Karzanov's maximum network flow 

algorithm will correctly compute, in time 0(|V|-|E| + |V|-|V| 2 ) - 0(|V| 3 ), a 

maximum flow function on any network given as input. 
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2.G Further Improvements 

Ue saw in the last section that any 0(f) algorithm for computing a 
maximal flow function on a layered network could be used to develop an 
0(|V|-|E| + |V|-r) maximum network flow algorithm. (Simply modify 
Karzanov's max-flow algorithm by replacing his "push and balance" procedure 
with the new maximal flow procedure.) Notice, however, that any max-flow 
algorithm developed through this technique can be no faster than 
0(|V|-|E|). Despite this fact, each of the four mo9t recent maximum 
network flow algorithms have been based on developing new maximal flow 
procedures. Ue shall now briefly describe each of these new maximal flow 
procedures. 

In 197G B.V. Cherkasky [18] showed how Karzanov's 0(|V| 2 ) push and 
balance routine could be modified to run in time ( | V | • | E | 1 ' 2 ) . 
Cherkasky* s procedure partitions the layered network into blocks of 
consecutive layers called super layers. It then applies Karzanov's push and 
balance techniques to these super layers. Within the super layers, however. 
Dime's flow augmentation techniques are used. The result i9 an 
asymptotically faster but very complex maximal flow procedure. 

A little over a year later, Z. Galil [10] improved Cherkasky' s 

routine to obtain a maximal flow procedure with time complexity 

2/3 2/3 
0(|V| -|E| ). Gal it's procedure differs from Cherkasky's procedure in 

the techniques used within the super layers. Galil's routine maintains a 

special data structure containing information about the current 

"usefulness" of chains within the layered network. This data structure is 
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used to expedite the push of flow through edges within the superlayers. 
Like Cherkasky's routine, Galil's procedure is very complex and requires a 
great deal of overhead. 

n 

In 1978 a very simple 0(|Vp) maximal flow procedure was developed 
by three Indians named V.M. Malhotra, tt. Pramodh Kumar, and S.N. 
Maheshwari [161. Their procedure is based on successively augmenting an 
existing legal flow function on the layered network. The procedure begins 
by determining the maximum amount of flow that can be pushed through each 
vertex v in the layered network. This value is called the f low potential 
/» f (v) of the vertex v. Each flow augmentation is then performed in three 
steps. First a vertex with minimum non-zero flow potential over all 
vertices in the layered network is selected as the reference vertex, r. 
Next, /b f (r) units of flow are pushed from r to t and from s to r. The 
pushing is performed essentially as outlined in Karzanov's algorithm. 
Finally, the procedure updates the flow potential of each vertex through 
which flow has been pushed, closing all edges which become either saturated 
or unreachable from s or t. Although this procedure is not an asymptotic 
improvement over Karzanov's push and balance routine, it is extremely 
simple and results in perhaps the best max-flow algorithm for use on dense 
networks ( |E| ~|V| 2 ) . 

Finally, Gal i I and A. Naamad 112] have recently developed an 

9 

( I E I - 1 og I V | ) maximal flow procedure. This procedure is similar to 
Galil's 0(|V| 2/3 -|E| 2/3 ) maximal flow procedure in that a data structure is 
maintained for processing chains within the layered network. The new 
algorithm, however, does not partition the layered network into 
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super layers. 

This concludes our discussion of algorithms for computing a maximum 
flow function on a network. We will now turn our attention to the question 
of lower bounds on the computational complexity of the maximum network flow 
problem. 



52 
CHAPTER 3 - LOWER BOUNDS ON MAX-FLOW 

3.1 Introduction 

Chapter 2 dealt with the establishment of upper bounds on the 
computational complexity of the maximum network flow problem. He traced 
the development of max-flow algorithms from the original labeling algorithm 
of Ford and Fulkerson [9], through the recent 0( |V| -|E| -log 2 |V | ) algorithm 
of Gal i 1 and Naamad [12]. The very fact that the search for new max-flow 
algorithms has been so fruitful leads us to now ask the question, "Can ue 
do better?". Can we develop a max-flow algorithm which is asymptotically 
faster than 0( |V| -|E|-log 2 |V| }? Gait I has shown [11,121 that any algorithm 
which uses Dinic's technique of dividing the problem into phases (as do all 
the known algorithms developed since Dinic's algorithm) must have time 
complexity at least 0(|V|-|E|). Further, he has conjectured an fl(|V|-|E|) 
lower bound on the computational complexity of the maximum network flow 
problem. At this time, however, there is no known non-linear lower bound 
on max-flow. 

The determination of lower bounds on the computational complexity 
of a problem is generally much more difficult than the establishment of 
upper bounds on the problem. In the latter case we can simply demonstrate 
an algorithm for solving the problem within the specified running time. In 
the former case, however, we must prove that any algorithm for solving the 
problem must require at least the specified running time, regardless of how 
clever the algorithm. In order to somewhat simplify the lower bound 
problem, many authors have chosen to work with restricted models of 
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computation. One such model which has received considerable attention in 
the recent literature is the linear decision tree model [6, 18,21,22] . This 
model tends to underestimate total time complexity but nevertheless enables 
us to study non-trivial lower bounds. In this chapter we shall investigate 
one particular approach to establishing non-linear lower bounds on the 
computational complexity of the maximum network flow problem relative to 
the linear decision tree model of computation. The technique we shall deal 
with is the polyhedral technique developed by A.C. Yao, D.fl. Avis and R.L. 
Rivest [21,221. 

3.2 The Node I of Computation 

The linear decision tree model of computation is based on the 
notion of a linear decision tree algorithm. A linear decision tree 
algorithm, operating on input (x^,...,* ), is simply a finite ternary tree 
with each internal node representing a test of the form "£ a j' x j :z " and 
each leaf containing a possible output. Given any input, the algorithm 
begins at the root and proceeds by moving down the tree until a leaf is 
reached. At each internal node the algorithm performs the specified test 
and then branches according to the result of this test (<,-,or >) . Once a 
leaf is reached, the information contained in that leaf is output as the 
result of the computation and the algorithm halts. The time complexity of 
any such algorithm is simply defined to be the height of the corresponding 
tree. 

The computational complexity of any problem relative to the linear 
decision tree model of computation can now be defined as the minimum height 
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over ail decision trees which solve the problem. It should be clear that 
the linear decision tree model measures time complexity solely in terms of 
the number of comparisons and branchings required. Thus the model tends to 
underestimate the total time complexity of a problem. Despite this fact, 
the linear decision tree model has proven useful in establishing several 
non-trivial lower bounds. 

3.3 Polyhedral Decision Problems 

Let P - {xelR n | & t {x)<d for each i-1,2 m} be a set of points in 

IR n , where 3»{xj, . . . ,x ) , m is an integer and 

n 

V*> ■ SXjj-Xj 

i-i 

for real numbers Xjj. The set P is said to be a polyhedron in R n . If d 
is the dimension of the smallest subspace of IR n containing P, then P is 
also said to be a polyhedron of dimension d. (Notice that we are 
restricting our attention to homogeneous polyhedra, i.e. cones.) On each 
subset H of the set {l,2,...,m}, we define the set of points F^P)cP as 
fol lows: 

F H (P) = {xelR n | l-(2)<e for each ieh\ 

ij (x)=0 for each i*H}. 
The set Fj_j (P) is called a face of the polyhedron P. I f s is the dimension 
of the smallest subspace of R n containing Fu(P), then F^(P) is said to be 
a face of dimension s. (The empty face has dimension -1 by convention.) Ue 
shall let F S (P) denote the set of all faces of dimension s of P. Notice 
that every point in P lies on some face of P and that the intersection of 
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any two faces of P is empty. 

The polyhedral dec i s i on problem B(P) can now be defined as the 
problem of determining whether an input point xeR 1 " 1 lies in the polyhedron 
P (i.e. Given any input xelR n , is xeP?). A linear decision tree algorithm 
for solving this problem will be a decision tree which contains a "yes" or 
"no" decision at every leaf. The computational complexity of B(P) relative 
to the linear decision tree model of computation will be denoted by L(P). 

In 1977 A.C. Yao, O.H. Avis and R.L. Rivest [21] proved the 
following fundamental theorem relating the complexity of the polyhedral 
decision problem B(P) to the facial character of the polyhedron P. 

Theorem 3.1 : 

Let P - {x | i.(x)<0 for each i-1,2 m} be a polyhedron in R n . 

Then for each s, 

L(P) > 1/2 log |F(P)|. 

Proof ; 

The proof of Theorem 3.1 can be found in C21J. 

This theorem states that fldog F s ) linear comparisons are necessary to 
determine if a point lies in a polyhedron composed of F_ s-dimensional 
faces. As a result, we can determine lower bounds on the computational 
complexity of any polyhedral decision problem relative to the linear 
decision tree model of computation by simply examining the facial structure 
of the polyhedron. 
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3.4 Applications to Max-Flow 

In this section we shall consider a straightforward application of 
the concepts presented in the previous section to the problem of 
establishing a non-trivial lower bound on max-flow. Ue shall first 
introduce the class of polyhedral decision problems {B(P n ) I n>2}, which is 
very closely related to the maximum network flow problem. Ue shall then 

9 

formalize this relationship by showing that L(P n )-n -2 is in fact a lower 
bound on L p , where l_ n is the linear decision tree complexity of the maximum 
network flow problem for a complete network on n vertices. Thus any lower 
bound on L(P n ) will also yield a lower bound on L . Finally we shall prove 
the fol lowing three results concerning the facial structure of the 
polyhedra P n (n>2) : 

1) There exists a positive constant c such that 
|F 1 (P n )| > c-(n-2H, for all n>2. 

2) There exists a positive constant c' such that 
|F 1 (P n )| < c'-(n-2H, for all n*2. 

3) There exists a positive constant c" such that 
|F g (P n )| < 2 (c "' n , , for all s and for all n>2. 

Based on these results, we can then conclude that Theorem 3.1 cannot be 

directly applied to the class of problems {B(P n ) | n*2} to obtain a non- 
trivial lower bound on max-flow. 

In order to formally define the class of polyhedral decision 

problems {B(P n ) | n>2}, we must introduce some new notation. Let G R denote 

the complete directed graph on ni2 vertices in which the vertex v^ is 

specified as s and the vertex v n is specified as t. Ue can represent the 
sets of vertices and edges of G n as follows: 

V - {vi.v2.v3 v n) 
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and 
E = { ejj | l<i,j<n). 
Notice that |V|=n, |E|=n and we are defining e : - to be the edge in G p 
directed from vertex v ( to vertex v.. It should be clear that any capacity 
function defined on the set E wi 1 1 give rise to a complete flow network on 
n vertices. Further, notice that the definition of a cut in a network can 
be directly applied to the graph G n . Ue now let {X^.Xo, . • . ,)(2 n - 2 } represent 
the set of all cuts in G . 

For our purposes, a vector y*eR n will be represented as 

-» n 2 +l 

y= ( yil yi n y nl•••••ynn' y n2+l , • The polyhedron P R in R n +1 can now 

be defined as follows: 

P n = {yelR n +1 | yjj >8 for all l<i,jsn, y n2+1 *0, ^(yJ^B for all l<k<2 n ~ 2 }, 
where 

V3> - (S {uij I «ij^tx k ;>< k >0-yn2 + i- 

The polyhedral decision problem B(P n ) is to determine whether an input 
point y*elR n belongs to the polyhedron P R . If we think of the set 

{y^l y nn5 as defin ' n g a capacity function on the graph G n (i.e. 

c(ej.)-yj.), then it should be clear that a point ycR n will belong to 
the polyhedron P n if and only if the following two relations hold: 

1) y ( ->8, for all l<i,j<n (i.e. c (e j s J "U j : defines 
a legal capacity function on G n ) . 

2) ^yr^+i-CMmin). where X mjn is any minimum 
capacity cut on the network defined by G 
and c(ej :)=yj 

By Theorem 1.3, however, we have that C(X • ) is equal to the value of any 
maximum flow function f |nax on the network defined by G p and c (e j j) «y j j . 
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Therefore, it should also be clear that a point yeR n +1 will belong to the 
polyhedron P n if and only if c (e • .) =y j- . defines a legal capacity function 
on G n and there exists some legal flow function f on the network defined by 
G n and c(e j j)=y. | j, such that v(f)=y n2+1 (i.e. 8<y n2+1 <v ( f mg>< ) ) . The 
following lemma relates the linear decision tree complexity of B(P n ) to the 
linear decision tree complexity of the max-flow problem on a complete 
network of n vertices. 



Lemma 3. 1: 



L n > L(P n )-n 2 -2. 
n n 



Before presenting the proof of Lemma 3.1, we should consider the 
structure of any linear decision tree algorithm which computes a maximum 
flow function on a complete network on n vertices. Such an algorithm will 
be a ternary tree operating on input (y*|,...,y ), where the input defines 
a capacity function c(ejj)«y.. on the graph G n . Each leaf in the tree will 
contain a set of n linear functions {g:. | l<i,j<n} defined on the n 
input variables. For any input, the algorithm will begin at the root and 
proceed by moving down the tree until a leaf is reached. Once a leaf is 
reached, a maximum flow function f on the network defined by the graph 

G n and the input (y^ y nn ) , will be given by f nax ^}-^ ' 

^i i * y ll' * ' ' *^nn^ ' ^ e snou 'd note that on every network there will exist a 
maximum flow function that can be completely defined by a set of linear 
combinations of the edge capacities on the network. 
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Proof of Lemma 3.1 : 

Let T be any optimal linear decision tree algorithm for computing a 
maximum flow function on a complete network on n vertices. Clearly the 
height of T must be L R . Further, we can obtain a linear decision tree T' 
for the problem B(P n ) by modifying T as follows. Place the root of T belou 
a new sequence of n 2 distinct tests of the form "Is yj j>8?" such that the 
root of T is reached if and only if all of these new tests produce a "yes" 
answer. Then replace each leaf in T with a new test of the form "Is 
y n2+1 >8?" followed by a new test of the form "Is v(f ma>< )-y n2+1 >8?" (y p2+1 
being a new input). The new tree T" should be constructed in such a way 
that if any of the newly added tests produce a "no" answer, then a leaf 
containing a "no" decision is reached. Otherwise a leaf containing a "yes" 
decision is reached. Since the value v(f mg>< ) is simply a sum of the gj . 
available at each leaf in T, it should be clear that T' is in fact a linear 
decision tree algorithm for solving the problem B(P n ). The height of T', 

9 

however, is L n +n^+2 and thus we have the relation: 

L(P n ) < L n+ n 2 +2 



L n > L(P n )-n 2 -2. 



□ 



The problem of establishing a non-trivial lower bound on the linear 
decision tree complexity of the maximum network flow problem has now been 
reduced to the problem of establishing a non-trivial lower bound on L(P p ). 
for each n>2. By Theorem 3.1, however, we can establish lower bounds on 
L(P n ), for each n>2, by simply examining the facial structure of the 
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polyhedra P n . The remainder of this section Mill be devoted to proving 
several lemmas concerning the number of faces composing each of these 
polyhedra. Ue wilt essentially show that each polyhedron is composed of 
relatively few faces and thus Theorem 3.1 can be of no use in establishing 
a non- linear lower bound (if one exists) on any L(P ). 

Lemma 3.2 : 

There exists a positive constant c such that |P^P n M * c-(n-2)!, 
for a I I n>2. 

Proof ; 

Consider any n>2 and let p' be an s-t chain in the graph G n . Now 
consider the point y'elR n +1 with the properties: 

1) yf:=l if ejjEp', for all Isi.jSn 

2) y( j= if ejj*p', for all l<i,j<n 

3) y n2+1 -i. 

Clearly the point y' belongs to the polyhedron P and so there must exist 
some face of P R containing y' (since every point in a polyhedron lies on 
some face of that polyhedron). Let Frt'(P p ) denote the face of P n which 
contains the point y'. Therefore, F^(P p ) is the set of all points y"eP n 
which satisfy: 
<3.1) y(j=0 iff y[j=0, for all lsi,j<n 

(3.2) y^ +1 >0 

(3.3) i^jJ'-'J-e iff i k (y')-0, for all lsks2 n_2 . 

From (3.1) and (3.2) it should be clear that every point y"eF^(P n ) will 
be non-zero in exactly the same co-ordinates. Further, applying this 
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observation to (3.3) and the fact that p' is an s-t chain in G n , we must 
have that M\'\=\Arf.+\ f° r eacn non-zero yj^ (l<i,j<n). Thus it must be the 
case that every point in the set Fj->/(P n ) is simply a scalar multiple of the 
point y'. It should now be clear that the smallest subspace of R 
containing the set Ftf'(P n ) has dimension 1. Thus we have that 
F^-'(P n ) eFj (P n ) . Finally, notice that each distinct s-t chain in G n will 
give rise to a distinct 1-dimensional face of P n> Therefore, if Z denotes 
the number of distinct s-t chains in the graph G n then we must have: 

(3.4) l F l (P n } ' * Z * 

Since G p is a complete graph on n vertices, however, the value Z can be 

expressed by the following formula: 

n-1 

Z - 2 (n-21 !/(!-»! (8! - 1) 

i-1 

n-2 

- (n-2)! -2 1/H 
i-8 

> (n-2)! 

Combining this expression with (3.4) we obtain: 

IF^PJI > (n-2)! 



Lemma 3.3 : 

There exists a positive constant c' such that |Fj(P n )| £ c'-(n-2)!, 
for a I I n>2. 

Proof : 

Consider any n>2 and let F{(P n ) denote the set of all 1-dimensional 
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faces of the polyhedron P n which arise from distinct s-t chains in the 
graph G n , as outlined in the proof of Lemma 3.2. Now let Q denote the set 
of all points $'eP n which satisfy any one of the following three 
condi t i ons: 
1) y'=8 

2 ^ ^n2+l = ®' yfn = ® anc * 3' ^ as exac ^'y one non-zero co-ordinate 
3) y' belongs to a face in the set Ff(P ). 

First notice that each point in Q which satisfies condition (2) belongs to 

2 
one of exactly n -1 distinct (trivial) 1-dimensional faces of P n . Further, 

9 

notice that none of these n -1 faces belongs to the set Ff(P n ). Thus each 
point in Q, except the point 8, belongs to one of exactly |F£ (P p > |+n^-l 
distinct 1-dimensional faces of P p . Next observe that every point in the 
polyhedron P n can be expressed as a convex combination of points in Q 
(consider expressing a network as a sum of distinct s-t chains and isolated 
edges). Therefore, every point on a 1-dimensional face of P R can be 
expressed as a convex combination of points in Q. Thus we must have that 
the set Q contains at least one point belonging to each 1-dimensional face 



of P and hence: 



f p l {P n ) > * l F l {P n , l +n 2 - 1 
n-2 

> ((n-2)! X 1/i! ) + n 2 - 1 
i-8 

£ 3-(n-2)! + n 2 - 1 

s c'-(n-2)! (for some positive constant 

c", independent of n) 



D 
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LemmS 3.4: 



There exists a positive constant c" such that |F_(P_)| < 

9 n 



2 VC n ', for all s and for all n>2. 



Proof : 

,2 



Consider any n>2 and let each vector weR be represented as 

-♦ ■ — - 7 2 

w=(w^, . . . ,w nn ,ui^,. .. ,w pn ) . A polyhedron P^ in fr n can now be defined 



relative to the graph G as follows: 



,2 



P^ = {welR 2rr | Wj.>8 and Wj.>0 and w-.iw. . for all l<i,j<n, 

£'(w)>0 and jE'(i3)<0 for all 2<q<n-l, 
J£'(iJ)*8 for all l<k<2 n ~ 2 }, 
where 

^q<"> - Z {"jj I e ije ln(v q )} - £ {«,. | ejj 60ut(v q )} 
and 

^'(w) - £ { Wij | e ije (K k; R k )} - 2 {u. (j | e|j eIn(t)J. 
If we think of the set {wi , , . . . f w } as defining a capacity function on the 
graph G p (i.e. c(e ( .)=Wj.) and the set {w^, . . . »w" nn } as defining a flow 

function on the network defined by G n and c(ej.)-w.. (i.e. f (ej )-Wj •) , 

7 2 
then it should be clear that a point weJr will belong to the polyhedron 

P^ if and only if the following two conditions hold: 

1) c(ej.)=Wj- defines a legal capacity function 
on the graph G n (i.e. c(e)>0 for each eeE). 

2) f(e ; ;)=w ; . defines a legal flow function on 
the network defined by G and c(e- :)-w. .. 

Notice that the constraints ££'(3)20 are all redundant since any legal flow 
function f on a network Tl must always satisfy the relation v(f )<C(X m j n ) , 
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where X mm is any minimum capacity cut on H. It should now be clear that 
we can determine whether a point weJr n belongs to the polyhedron P^ by 
simply testing to see if d satisfies each of the first 3n 2 +2n-4 constraints 
defining P^. Thus there exists a straightforward linear decision tree 
algorithm of height 3n^+2n-4 for solving the problem B(P^) and hence we 



have the relation: 



L(P n ) < 3n 2 +2n-4 < 4n 2 . 



Further, combining this relation with Theorem 3.1 we obtain the result: 

4n 2 > L(P n ) i 111 log |F g (P n )|, for all s 

(3.5) => 2 8n > |F s (P n )|, for all s. 

If we now let F(P') -■ U FJP') denote the set of all faces of the 

" s 

polyhedron P' then by (3.5) we have: 



|F(P n )| < £2 



8n 2 



-l*s*2n 2 



=» |F(P^)| < (2n 2 +2)-2 8n 

(3.G) => |F(P n )| < 2 {c "" n2) 

(for some positive constant c", independent of n) . 

Turning our attention back to the polyhedron P R , let F(P p ) - 
U F s (P n ^ denote the set of all faces of the polyhedron P n . The remainder 
of this proof will essentially consist of constructing a 1-1 mapping \ff 
from the elements of F(P n ) into the elements of F(P^). He will then have 
that |F(P n )| < |F(P^)|. 

Let F H (P n )eF(P n ) be any face of the polyhedron P n . If F H (P n )-0 
(i.e. F^j(P n ) is the empty face) then define ^(F^(P n ) ) «=$. If, however, 
F H (P n )*^ then let y"elR n +1 be some point in the set F H (P n ). Since 
y'eF H (P n )cP n , we must have that c(ej.)-y^ defines a legal capacity 
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function on the graph G n and that there exists some legal flow function f' 

on the network defined by G n and cle^-Ny^'-, such that v ( f ' ) -y^i • Now 

7 2 
consider the point w'eR such that: 

wfj-yfj, for all lsi, j<n 

and 

w[.-f'(ej.), for all l<i , j<n. 

Clearly w'eP^ and thus there must exist some face of P' containing w". 

Call this face F^/IP^) and define ^(F H (P n ) NF^/ (P^) . It should now be 

clear that ^ maps each face in F(P n ) into a face in F(P^). Thus it 

remains to be shown only that ^ is a 1-1 mapping. Recall, however, that 

each face of a polyhedron is uniquely determined by the set of constraints 

its elements satisfy by equality. Further, notice that the points w' and 

y' have the following relationship: 

1) w(.=e iff y(:-8, for all l<i,j<n 

2) ^'(w')-8 iff & k $')*d, for all lsks2 n " 2 

3) w(j=0 for all {i j | e, ,eln(t)} iff yn2 + i= - 

Thus for each distinct face F H (P p ) of the polyhedron P n , ^(F^|(P n )) will be 

a distinct face of the polyhedron P^ and so ^ is in fact a 1-1 mapping. 

Ue can now conclude that |F(P n )| < |F(P p )| and thus by (3.G) we have: 

(3.7) |F(P n ) | < 2 (c "* n2) . 

Finally it should be clear that |F g (P p )| i |F(P n )|, for all s, and so by 

(3.7): 

|F g (P n )| i 2 (c "' n2) , for all s. 
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3.5 Conclusions 

In this chapter we presented the linear decision tree model of 

computation [6,18,21,22], the notion of a polyhedron and a polyhedral 

decision problem [21,22], and the class of polyhedral decision problems 

(B(P n ) | n>2} which most naturally arises when considering the maximum 

network flow problem. Ue then showed that the problem of establishing a 

non- linear lower bound on the linear decision tree complexity of max- flow 

can be reduced to the problem of establishing a non-linear lower bound on 

l_(P n ), for each n>2. Next we demonstrated matching upper and lower bounds 

on the number of faces of dimension 1 composing each of the polyhedra P p 

2 
(n>2). Finally, we established a 2° ' upper bound on the number of 

faces, of any dimension, composing each of the polyhedra P n . Based on our 

results, we can now conclude that Theorem 3.1 can be of no use in 

establishing a non-linear lower bound on max-flow, through the class of 

polyhedral decision problems (B(P n ) | n>2}. It remains an open question, 

however, whether or not the polyhedral technique can in general be useful 

for establishing non-linear lower bounds on max-flow. There may, for 

example, exist some more complex class of polyhedral decision problems that 

can be reduced to max-flow. 
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